LAPLACIANS ON PERIODIC DISCRETE GRAPHS 



ANDREY BADANIN, EVGENY KOROTYAEV, AND NATALIA SABUROVA 



Abstract. We consider Laplacians on Z 2 -periodic discrete graphs. The following results 
are obtained: 1) The Floquet-Bloch decomposition is constructed and basic properties are 
derived. 2) The estimates of the Lebesgue measure of the spectrum in terms of geometric 
parameters of the graph are obtained. 3) The spectrum of the Laplacian is described, when 
the so-called fundamental graph consists of one or two vertices and any number of edges. 4) 
We consider the hexagonal lattice perturbed by adding one edge to the fundamental graph. 
There exist two cases: a) if the perturbed hexagonal lattice is bipartite, then the spectrum 
of the perturbed Laplacian coincides with the spectrum [—1,1] for the unperturbed case, b) 
if the perturbed hexagonal lattice is not bipartite, then there is a gap in the spectrum of the 
perturbed Laplacian. Moreover, some deeper results are obtained for the perturbation of the 
square lattice. 



1. Introduction and main results 

1.1. Introduction. We discuss the spectrum of Laplacians on periodic discrete graphs. Lapla- 
cians on periodic graphs are of interest due to their applications to problems of physics and 
chemistry. They are used to describe and to study properties of different periodic media, 
including nanomedia, see |NG04| . [Ha85] . |SDD98j . 

There are a lot of papers, and even books, on the spectrum of the discrete Laplacian on an 
infinite graph. One of the main problems is to describe the spectral properties of the Laplacian 
in terms of geometric parameters of the graph. Useful sources concerning operators acting 
on infi nite gra p hs are t he book s [Ba98]. jBo^8] . |BH12| . |GDS95] . |GDGT88] . [P12] and the 
papers |HN09j . |Me94j . |MW89j . [MRA07] . See also the references therein. 

There are papers about the spectrum of discrete Laplacians on periodic graphs. Higuchi 
and Shirai [HS04] (see also [RR 07] ) obtain the decomposition of the Laplacian into a constant 
fiber direct integral. Higuchi and Nomura [HN09J prove that the spectrum of the Laplacian 
consists of an absolutely continuous part and a finite number of flat bands (i.e., eigenvalues 
with infinite multiplicity). The absolutely continuous spectrum consists of a finite number of 
intervals (spectral bands) separated by gaps. Moreover, they also show that for each flat band 
there exists a finitely supported eigenfunction. 

There are results about spectral properties of the discrete Laplace and Schrodinger opera- 
tors on specific periodic graphs. The hexagonal lattice can be viewed as a discrete model of 
graphene, which is two-dimensional single- layered carbon sheet with honeycomb structure. It 
was recently discovered by Geim and Novoselov |NG04] . Graphene is a very hot subject in 
physics, where tight binding models for the Schrodinger operator are standard and give inter- 
esting band structures. Ando |A12j considers the spectral theory for the discrete Schrodinger 
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operators with finitely supported potentials on the hexagonal lattice and their inverse scat- 
tering problem. Korotyaev and Kutsenko |KK10] - |KK10b] study the spectra of the discrete 
Schrodinger operators on graphene nano-tubes and nano-ribbons in external fields. See more 
about graphene in Section 5. Schrodinger operators with decreasing potentials on the lattice 
Z d are considered by Boutet de Monvel-Sahbani [BS99], Isozaki-Korotyaev |IK12] . Rosenblum- 
Solomjak |RoS09] . Gieseker-Knorrer-Trubowitz |GKT93j consider Schrodinger operators with 
periodic potentials on the lattice Z 2 , the simplest example of Z 2 -periodic graphs. They study 
its Bloch variety and its integrated density of states. 
We describe our main goals of the paper: 

1) to construct the Floquet theory for Laplacians on periodic graphs. 

2) to estimate the Lebesgue measure of the spectrum of the Laplacian in terms of geometric 
parameters of the graph. 

3) to describe the spectrum of the Laplacian under the perturbation of the square lattice 
and the hexagonal lattice. 

4) It is known |C97j that the investigation of the spectrum of the Laplacian on an equilateral 
quantum graph (i.e., a graph consisting of identical segments) can be reduced to the study of 
the spectrum of the discrete Laplacian. Thus, we have to do the needed spectral analysis of 
the discrete Laplacian given in this paper, in order to describe spectral properties (including 
the Bethe-Sommerfeld conjecture) of the Laplacians on quantum graphs [KSJ. 

It should be noted that the results obtained in this work can be generalized to the case of 
Z d , d ^ 3, periodic graphs |KS1] . 

1.2. The definitions of periodic graphs and fundamental graphs. Let T = (V,£) be a 
connected graph, possibly having loops and multiple edges, where V is the set of its vertices 
and £ is the set of its unoriented edges. The graphs under consideration are embedded into 
M 2 . Considering each edge in £ to have two orientations, we can introduce the set A of all 
oriented edges. The inverse edge of e G A is denoted by e. The oriented edge starting at 
u G V and ending at v G V will be denoted as the ordered pair (u,v). Vertices u,v G V will 
be called adjacent and denoted by u ~ v , if (u,v) G A. We define the degree x v of the vertex 
v G V as the number of all oriented edges from A starting at v. Below we consider Z 2 -periodic 
graphs T, satisfying the following conditions: 

1) the number of vertices from V in any bounded domain C M 2 is finite; 

2) the degree of each vertex is finite; 

3) there exists a basis ai,a 2 in M. 2 such that T is invariant under translations through the 
vectors a\ and a 2 : 

r + a\ = T, T + a 2 = T. 

The vectors a±,a 2 are called the periods of'T. 

In the plane 1R 2 we consider a coordinate system with the origin at some point O. The 
coordinate axes of this system are directed along the vectors a\ and a 2 . Below the coordinates 
of all vertices of V will be expressed in this coordinate system. Then it follows from the 
definition of Z 2 -periodic graph that V is invariant under translations through any integer 
vector: 

r + p = r, v^gz 2 . 

Examples of periodic graphs are shown in Figures H]a, [7]a. 
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We define the fundamental graph T of the periodic graph T as a graph on the surface M 2 /Z 2 
by 

T = T/Z 2 C T 2 , where T 2 = R 2 /Z 2 . (1.1) 

The vertex set Vb, the set So of unoriented edges and the set Aq of oriented edges of r are 
finite (see Proposition 16.21) . 

We introduce an edge "index", which is important in the study of the Laplace operator. We 
identify the vertices of the fundamental graph T = (Vq, £q) with the vertices of the periodic 
graph T = (V,£) in the set [0, l) 2 . Then for any v G V the following unique representation 
holds true: 

v=[v] + v, [v] G Z 2 , veV C [0, l) 2 . (1.2) 

In other words, each vertex v can be represented uniquely as the sum of an integer part 
[v] G Z 2 and a fractional part v that is a vertex of the fundamental graph IV For any oriented 
edge e = (u,v) E A we define the edge "index" r(e) as the integer vector 

r(e) = [v] - [u] G Z 2 , (1.3) 

where 

u = [u] + u, v = [v ] + v, [u], [v] G Z 2 , U, V G Vq. 

If e = (u, v) is an oriented edge of the graph T, then by the definition of the fundamental 
graph there is an oriented edge e = (u, v ) on r . For an edge e G we define the edge index 
r(e) by 

r(e) = r(e). (1.4) 
In other words, indices of periodic graph edges are inherited by edges of the fundamental 
graph. The edge indices, generally speaking, depend on the choice of the coordinate origin 
O. But in a fixed coordinate system the index of the fundamental graph edge is uniquely 
determined by (11.4)) . since due to Proposition I6.3l iii we have 

r(e + p) = r(e), V (e,p) G -4 x Z 2 . 

Edges with nonzero indices will be called bridges. They are important, when we describe the 
spectrum of the Laplacian. The bridges provide the connectivity of the periodic graph and 
the removal of all bridges disconnects the graph into infinitely many connected components. 

1.3. Laplace operators on graphs. Let £ 2 (V) be the Hilbert space of all square summable 
functions / : V — > C, equipped with the norm 



2 _ I f/„.\\2 

p-iy) — 



Recall that the degree h v of the vertex v G V is the number of oriented edges starting at the 
vertex v. We define the Laplacian (or the Laplace operator) A acting on the Hilbert space 
£ 2 (V)bj 

(a/)m = -L= J2 -j^f^ v/G£ 2 (n (i-5) 



v (v,u)eA v 

Note that sometimes 1 — A ^ is also called the Laplace operator. 

A graph is called bipartite if its vertex set is divided into two disjoint sets (called parts 
of the graph) such that each edge connects vertices from distinct sets (see p. 105 in [Or62j). 
Examples of bipartite graphs are the square lattice (Figure Wp) and the hexagonal lattice 
(Figure Sfr). The triangular lattice (Figure [9|)) is non-bipartite. Note that for a bipartite 
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periodic graph there exists a bipartite fundamental graph (see Lemma [6.4l ii). but not every 
fundamental graph is bipartite. Indeed, the square lattice (Figure 0a) is bipartite, but its 
fundamental graph shown in Figure ^7$ is non-bipartite. 

We recall well-known properties of the Laplacian A, which hold true for finite and infinite 
graphs (in particular, for periodic graphs) (see |Ch97] . |Me94] . |M92] , |MW89] ) : 

Main properties of the Laplacian: 

1) The operator A is self-adjoint and bounded. 

2) The spectrum a (A) is contained in [—1, 1]. 

3) The point 1 belongs to the spectrum cr(A). 

4) The graph is bipartite <=> the point — 1 G er(A) ^ the spectrum cr(A) is symmetric with 
respect to the point zero. 

Denote by Vi, . . . , v v the vertices of Vb, where v < oo is a number of vertices of Vo- Thus we 
have 

V = {v j :jeH v }, N„ = {l,...,i/}. 
We present our first results about the decomposition of the Laplacian A into a direct integral. 

Theorem 1.1. i) The operator A acting on £ 2 (V) has the decomposition into a constant fiber 
direct integral 



(27F) J * (1.6) 



for some unitary operator U . Here the fiber space £ 2 (Vo) = C u and A($) = {Ajk($)}, $ G T 2 , 
is a Floquet v x v matrix given by 

{ 7= E e'™'*, if (v 3 ,v k )eAo 
A jk (tf) = t ' e=( "^ t)e ^ , (1.7) 

I 0, if (vj,v k ) A 

where x,- is the degree ofvj and (• , •) denotes the standard inner product in IR 2 . In particular, 

Yl cos(r(e), if a loop (vj,Vj) G Ao 

0, if {vj,Vj) <£ Ao 

ii) Let A^(i?) be the Floquet matrix for A defined by fli. 7| ) in another coordinate system 
with the origin 0\. Then the matrices A^ifi) and A($) are unitarily equivalent for each 

deT 2 . 

Hi) The spectrum (j(A(i?)) is contained in [—1, 1] for each •§ G T 2 . 

iv) The entry A jfc (-) is constant iff there is no bridge (vj,v k ) G Ao- 

v) The point 1 is never a flat band of A. 

vi) The Floquet matrix A(-) has at least one non-constant entry Aj k (-) for some j ^ k. 

vii) A fundamental graph Tq is bipartite iff the spectrum cr(A(i?)) is symmetric with respect 
to the point zero for each $ G T 2 . 
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Figure 1. A periodic graph T; only edges of the fundamental graph T are 
shown; the bridges of T are marked by bold lines. 

Remark. 1) The representations (11.71) . fll.Sp are new. Recall that the existence of A($) was 
proved in |HS04] . 

2) The graph T with v = 5 on Figure Q] has the following bridges 

(vi,v 2 + ai), (vx,v 3 + a 2 ), + (^3, v A + a ± + a 2 ) 

because they connect vertices with different integer parts in the sense of the identity (11. 2p . 

3) There are always some bridges on a graph T, if T is connected. Otherwise, if there are no 
bridges on some graph T, then this graph is not connected and is a union of infinitely many 
connected components. In this case the spectrum of the Laplacian on T consists of only a 
finite number of eigenvalues with infinite multiplicity. Due to Theorem 11.41 i. it is impossible. 
The presence of bridges on the graph "enlarges" the spectrum of the Laplacian. 



1.4. Spectral properties of Laplacians. Theorem 11.11 and standard arguments (see The- 
orem XIII. 85 in [RS78]) describe the spectrum of the Laplacian. Each Floquet v x v matrix 
A($), G T 2 , has v eigenvalues \ n ($), n = I,..., v. They are real and locally analytic 
functions on the torus T 2 , since A($) is self-adjoint and analytic in $ on the torus T 2 . The 
parameter $ is called quasimomentum. If some A n (-) = C n = const on some set B C T 2 
of positive Lebesgue measure, then the operator A on V has the eigenvalue C n with infinite 
multiplicity. We will call C n a flat band. Due to |HN09j we have that A* is an eigenvalue of 
A iff A* is an eigenvalue of A($) for any d G T 2 . Thus, we can define the multiplicity of a flat 
band by: an eigenvalue A* of A has the multiplicity m iff A* = const is an eigenvalue of A(i?) 
for each d G T 2 with the multiplicity m (except maybe for a finite number of d G T 2 ). Thus 
if the operator A has r ^ flat bands, then we denote them by 

Hi = X v - j+1 (-&), j G N r , 

and they are labeled by 

(J>1 < f*2 < • • • ^ f*r, (1.9) 

counting multiplicities. Thus, all other eigenvalues A n (i?), n — 1, . . . , v — r are not constant. 
They can be enumerated in decreasing order (counting multiplicities) 

A I/ _ r (tf)<A„_ r _i(tf)<...<Ai(#), VtfGT 2 . (1.10) 
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Each A n ($) is a piecewise analytic function on the torus T 2 and defines a dispersion relation. 
Its plot is a dispersion curve. Define the spectral bands o~ n by 

cr„ = [AT, A+] = A„(T 2 ), where A~ = min X n 0&), A+ = max A n 0&), VnG N v . (1.11) 

For each n = l,...,i/-rwe have A~ < A+ and the spectral band a n is open (non-degenerate). 
For each n — v — r + 1, . . . ,v we have A^ = fi r and the spectral band a n is close (degenerate). 
If A^" +1 < A~ for some n G N„_ r _i, then the interval (A* +1 , A~) is called a gap. 
Thus the spectrum of the operator A on the periodic graph V has the form 

o-(A) = <7 ac (A)Uo7 6 (A), 

r: (i.i2) 

CT ac (A) = |J O n , 0-f b (A) = {[Mi, fi r }, 

71=1 

where cx ac (A) is the absolutely continuous spectrum and oy^A) is a set of all flat bands 
(eigenvalues with infinite multiplicity). We now describe precisely all bands for specific graphs. 

Theorem 1.2. i) Let all bridges of the fundamental graph T be loops, i.e., the indices of all 
edges (vj,Vk), 1 ^ j < k ^ u, ofT be zero. Then 

A+ = A n (0), VneN„. (1.13) 

Moreover, the eigenvalue A n (0) is a flat band of A iff A~ = A+. The multiplicity of the flat 
band \i = A^ of the Laplacian A is the multiplicity of A+ as the eigenvalue of the operator 
A(0). 

ii) Let, in addition, cos(r(e), i? ) = — 1 f or a ^ bridges e E £ and some i9 € T 2 . T/ien 

^n = [A;,A+], A; = A n (^ ), V n G Ny. (1.14) 

Lei T 6e a bipartite periodic graph, satisfying the condition i) (T is not bipartite, since 
there is a loop on Tq). Then A^ are the eigenvalues of the matrix ±A(0). 

Remark. 1) A+, n G N u , are the eigenvalues of the operator A(0) defined by (II. 7p that is 
the Laplacian on the fundamental graph To. The case of item ii) is similar to the case of 
iV-periodic Jacobi matrices on the lattice Z (and for Hill operators). The spectrum of these 
operators is absolutely continuous and is a union of spectral bands, separated by gaps. The 
endpoints of the bands are so-called 2iV-periodic eigenvalues. 

2) For some classes of graphs one can easily determine $o such that cos(r(e), i? ) = —1 for 
all bridges e G Aq. Let r(e) = (ri(e), r 2 (e)) G Z 2 be the index of a bridge e of r Q . Then 

{(ir, 0), if Ti(e) is odd for all bridges e E Ao 

(0, 7r), if r 2 (e) is odd for all bridges e G Ao . (1-15) 

(tc,tc), if Ti(e) + r 2 (e) is odd for all bridges eG^o 

3) The fundamental graph T on Figure [2]a has only 4 oriented bridges, which are loops with 
indices (0, ±1) and (±1, 0). All other edges of T have zero indices. Thus, 7~i(e) + r 2 (e) is odd 
for all bridges e G Ao and according to Theorem ll.2l i-ii and the identity (I1.15P the spectrum 
of the Laplacian on T, see Figured, has the form a (A) = Un=i a n, o~ n = [A n (7r, ir), A„(0)]. 

We estimate the spectrum a (A) of the Laplacian. 
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(a) 



(b) 




Figure 2. a) The fundamental graph T ; b) the periodic graph T. 
Theorem 1.3. i) The Lebesgue measure \a(A)\ of the spectrum of the Laplacian A satisfies 

bjk (1.16) 



u—r 



71=1 



where bjk is the number of bridges (vj,Vk) on the fundamental graph Tq, is the degree of 
Vj. Moreover, this estimate is sharp, i.e., estimates II 1.16]) become identities for some graphs. 

ii) For any v ^ 2 there exists a periodic graph such that |cr(A)| = In particular, 

|cr(A)| as v ->■ oo. 

Lei a fundamental graph T be bipartite. If v > 1 is odd, i/ien fi = is a flat band of A. 



Remark. 1) The estimate f II . 1 6 j) is effective for the case 

b jk 



< i. 



;i.i7) 



If the inequality (11.171) does not hold true, then (I1.16P yields only a simple estimate |cr(A)| ^ 2, 
which follows from the basic property er(A) C [—1, 1]. The inequality (11.17)) holds true when 
for each vertex v G Vq the ratio of the number of bridges starting at v to the degree of v is 
rather small. Increasing the degree of each vertex and fixing the bridges we can construct a 
graph such that the Lebesgue measure of the spectrum of the Laplacian is arbitrarily small. 

2) The total length of spectral bands depends essentially on the bridges on To- If we remove 
the coordinate system, then the number of bridges on the corresponding fundamental graph 
r is changed in general. 




Figure 3. The cycle that is the support of the eigenfunction with the eigen- 
value 0. The values of the eigenfunction in the vertices of the support are 0, ±1. 
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3) Recall that for each flat band there exists a finitely supported eigenfunction. It is well- 
known that for some specific case these eigenfunctions can easily be determined using some 
analysis of graph cycles. For example, if the graph has a cycle of length An, n G N, and the 
degree of each second vertices of the cycle is 2, then one can easily construct an eigenfunction 
supported on this cycle (see Figure [3]) with the eigenvalue 0. Theorems 11.21 and 1 1 . 3 1 iii give the 
simple sufficient conditions for the existence of flat bands, not based on an analysis of cycles. 
This is the only test for the existence of flat bands, that we know of. 

Now we describe possible positions of flat bands. 

Theorem 1.4. i) The first spectral band o\ = [A^,l] is open, i.e., Ajf < 1. Moreover, the 
number of flat bands r < v and the spectrum of the operator A has at most v — r — 1 gaps, 
ii) Let v ^ 2. Then 

1) There is a graph Y, such that the spectrum of the Laplacian on Y has exactly 2 open 
spectral bands and between them, in the gap, v — 2 flat bands, counting multiplicity. 

2) There is a graph Y, such that the spectrum of the Laplacian on Y has [|] different simple 
flat bands embedded in the absolutely continuous spectrum [—1,1]. 

3) There is a graph Y , such that a point — | is a flat band of the Laplacian on Y (the Kagome 
lattice, see subsection \5.5\) and lies at the endpoint of the absolutely continuous spectrum 
a ac (A) = [-i,l]. 

There is an open problem: does there exist a graph with any v ^ 2 vertices in the funda- 
mental graph such that the spectrum of the Laplacian on Y has only 1 spectral band and v — 1 
flat bands, counting multiplicity? 

1.5. Perturbations of the hexagonal lattice. Let G = (V,S) be the hexagonal lattice 
(FigureHJa). The periods of G are the vectors a\ = (3/2, y/3/2), = (0, v3) (the coordinates 
of ai, a2 are taken in the orthonormal basis ex, e^). The vertex set and the edge set are given 
by 

v = z 2 u(z 2 +(!,i)), 
£ = {(p,p+ (§,§)), (p,p+ Cp,p+ (§,-!)) v p ez 2 }. 

Recall that the coordinates of all vertices are taken in the basis ai, a<i- The fundamental 
graph Go consists of two vertices v±, t>2, multiple edges ei = e 2 = e 3 = (t>i,t>2) (Figure©) 
and ei = §2 = e3 with the indices 

r(ei)=r(ei) = (0,0), r(e 2 ) = -r(e 2 ) = (1, 0), r(e 3 ) = -r(e 3 ) = (0, 1). 
It is known that the spectrum of the Laplacian A on G has the form ex (A) = cr ac (A) = [—1,1]. 



V 2 




Figure 4. a) Graphene G; b) the fundamental graph Go of the graphene. 
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We describe the spectrum under the perturbation of the graphene G by adding an edge 
to its fundamental graph. 

Theorem 1.5 (Perturbations of the graphene). Let G' be a perturbed graph obtained from 
the hexagonal lattice G by adding one edge (i.e., two oriented edges e, e with the indices 
r(e) = — r(e) = (ti,T2)) to its fundamental graph Go- Then the spectrum cr(A') of the 
Laplacian A' on the perturbed graph G' satisfies: 

i) o-(A') = a ac (A'), i.e., the Laplacian on G' has no flat bands. 

ii) er(A') = [-1,1] iffG' is bipartite. 
Hi) Let G' be non-bipartite. Then 

tr(A') = [Aj,0]U[Ar,l], 0<Ar^^, z/ n - r 2 G 3Z, (1.18) 

( t(A / ) = [A 2 ,A+]U[0,1], -i^A+<0, ifn-T 2 ^3Z, (1.19) 

where 

-KA 2 ^-jj. (1.20) 
Moreover, the Lebesgue measure of the spectrum satisfies 

\a(A')\>-, if ri-r 2 G3Z; 

7 (1-21) 

|a(A')|^-, if n-r 2 ^3Z. 
5 

Remark. This theorem shows that if the perturbation keeps the bipartition of the graph, 
then the spectrum of the operator does not change. If the added edge breaks the bipartition 
of the graph, then there appears a gap in the spectrum and the lower point of the spectrum 
moves to the right. 

We present the plan of our paper. In section 2 we prove Theorem 11.11 about the decompo- 
sition of the Laplacian into a direct integral and prove some basic properties of the Floquet 
operator. In section 3 we estimate the Lebesgue measure of the spectrum in terms of geo- 
metric parameters of the graph in Theorem 11.31 In sections 4 and 5 we consider the periodic 
graphs with one, two and three vertices in the fundamental graph. In particularly, we describe 
the spectrum of the Laplacian on the square lattice and the hexagonal lattice, perturbed by 
adding edges to its fundamental graphs. In Appendix we recall some well-known properties 
of matrices and give some auxiliary statements, needed to prove our main results. 

2. Proof of Theorem 11.11 

In Appendix we recall some well-known properties of matrices, which will be used below. 
Also we prove that the fundamental graph of a periodic graph is finite and show that a 
bipartite periodic graph has a bipartite fundamental graph. Properties of an edge index, 
needed to prove our main results, are given in Appendix. 

Proof of Theorem fTTl i) Let = L 2 (T 2 , J^L , H) = J® 2 H A be a constant fiber direct 
integral, where "H = C u . It is the Hilbert space of square integrable functions / : T 2 — > C u , 
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equipped with the norm 



2 / II x/„QM|2 



Let y(V) be the set of all finitely supported functions / G £ 2 (V). Recall that v±, . . . ,v u are 
the vertices of the fundamental graph T . We identify the vertices of T with the vertices of 
the periodic graph T from the set [0, l) 2 . We define the operator U : £ 2 (V) — » M' by 

(Uf) n {#) = e- l{p ^f(v n +p), (#,7i) G T 2 x N„. (2.1) 

pez 2 

Using standard arguments (see pp. 290-291 in [RS78] ). we will show that U is well defined 
for / G y{y) and uniquely extendable to a unitary operator. 
For / G ^{V) the sum (12. ip is clearly convergent. For such functions / we compute 

f Hit 

\\ujt?= / \\{uf)mh> ' 



T2 (2tt)2 



/ E ( E e ~ i(M> /K ( E e * M> /K + ?)) 

n=l % e z 2 ' %ez 2 ' 

n=l dpZ 2 oeZ 2 ^ Jj2 



i(p-g,d) _ 



= E E I/k +^)l 2 = E i/(^)i 2 = ii/iij-oo- 

i=i pez 2 uev 

Here we have used the identity V = [v n + p : (n,p) £ N„ x Z 2 }. Then U is well defined and 
has a unique extension to an isometry. To see that U is onto Jf? we compute U*. We note 
that each vertex v G V can be written as 

v = v n + p for some (n,p) e N„ x Z 2 (2-2) 

and this representation is unique. Let g = (g n ($) : T 2 — > CJ ngN G J^ 5 . Then we define 

where p and n are defined by (12.21) . A direct computation shows that this is indeed the formula 
for the adjoint of U. Moreover, 



\U*9\\l {V) = EK^Mf = E E \(U*g)(v n +p)\ 2 = J2H f elM 9r 
nev n=l P ez 2 «=i pez 2 T2 



(■&) 



dd 2 



(2tt)' 



(27T) 

where we have used the Parseval relation for the Fourier series. 



\9\\je, 



LAPLACIANS ON PERIODIC DISCRETE GRAPHS 



11 



For / G y(V) and j G N u we obtain 



E^-LE E -^/(«*+P + r(e)) 



Xj A/Xt 

pez 2 v j t=i e =(« J ,« t )eAi v 



= y y — e'^-w y e -^ e ^> / + P + r (e)) 

fc=l e =(^,i) fe )e^o V J PGZ 2 

= E E ^=^^(f//) fc (^) = ^A, fc (^)(f//) fe (^), (2.3) 

where A jfc (0) are defined by (TTT7j) . The identity O yields (l/A/)(#) = A(tf)(l7/)(0), where 
A(#) = {A jfc ($)}. Thus, we obtain 

Using Proposition I6.3l i we can write the entries Ajj of the matrix A(i?) in the form (11. 8p . 
Thus, the statement i) has been proved. 

In order to prove ii) - vii), we discuss needed properties of the Floquet matrix A(0). 

Proposition 2.1. The matrix A(0) is the Laplacian on the fundamental graph T and has 
the form 

A(0) = {A jk (0)}, A jk (0) = V(j,fc)eC (2.4) 

where Xj k ^ 1 is the multiplicity of the edge (vj,v k ) G A$ and x^ = if (vj,v k ) ^ Ao- 
Moreover, they satisfy 

V 

Xj = ^2>c jk ^l, VjGN„. (2.5) 
k=i 

Proof. The identity (12. 4p for the matrix A(0) is obtained by a direct calculation of its entries 
using the formula (II. 7p . From the definition (II. 5p we deduce that the Laplacian on To has the 
form (12.41) . The definition of the degree Xj of the vertex Vj as the number of oriented edges 
starting at Vj gives (I2.5p . Since the graph T is connected, Xj ^ 1. ■ 

Proof of Theorem ll.ll ii — vii. ii) Recall that we identify the vertices V\, . . . ,v u of the 
fundamental graph r with the vertices of the periodic graph T from the set [0, l) 2 in the 
coordinate system with the origin O. Proposition I6.3l iv gives that for each (vj,v k ) G Ao 

r {1 \vj, v k ) = t(vj, v k ) + Pk - Pj, p k = [v k - b), pj = [vj - b), (2.6) 

where r^(vj, v k ) is the index of the edge (vj, v k ) in the coordinate system with the origin Oi, 
\. Using (12. 6p we rewrite the entries of A^(i5) defined by (ll.7p in the form 

A«(tf) = _L_ y e *<rW(e),*> = f y e ^ e ^ = e^-^A ik m. 

V J fc e=( Vj ,v k )eAo V 3 e=(*i, 
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We define the diagonal v x v matrix 

ZY(tf) = diag( e-*^'*), e-*<P">*> ) , Vi?gT 2 . 

A direct calculation yields 

U(#)A(i})ir 1 ($) = AW($), VtfeT 2 . 

Thus, for each $ 6 T 2 the matrices A^'(iJ) and A(i?) are unitarily equivalent. 

iii) Since A(0) is the Laplacian on T , ||A(0)|| = 1. From this fact and the formula (jl.7p it 
follows that the entries of each matrix A($), # G T 2 , satisfy 

|A ifc (#)| < Ajjfc(O) < 1, V (j, k) G N 2 . (2.7) 

Then, Proposition I6.1l i implies that the spectral radius p(A(i9)) satisfies 

p(A(0))<p(A(O)) = ||A(O)|| = l, 

which yields er(A(0)) C [-1, 1] for each G T 2 . 

iv) This statement follows immediately from (jl-7p and the definition of the bridge. 

v) We prove by the contradiction. Let the point 1 be an eigenvalue of the Laplacian on a 
graph T. Then there exists an eigenf unction ^ f G £ 2 (V) with eigenvalue 1 and with a finite 
support B C V (see Theorem 3.2 in |HN09j ). Let M = max = for some v G B. Thus, 

we have 

/(€) = (A/)(t;) ^= £ _L/(«)<-I= £ -±=/(€) = Mv^ = /(«)• 

We conclude that the inequality has to be an equality, and therefore 

/(«) /(«) v - 

a/ "y 

Repeating this argument until we reach a vertex from V \ B, we conclude that / = 0. We 
have a contradiction. 

vi) We prove by the contradiction. Let all entries Aj k (-), 1 ^ j ^ ^ v, be constant. Due 
to the self-adjointness of the matrix A($) all its entries are constants. Then the point 1 is an 
eigenvalue of A($) = A(0) for each G T 2 , i.e., the point 1 is an eigenvalue of A with infinite 
multiplicity. This contradicts item v). 

vii) Let To be a bipartite fundamental graph with the parts V\ = {t>i, . . . , v k } and Vi = 
{vk+i, ■ ■ ■ ,v u }. Since vertices from the same part of r are not adjacent, each matrix A(i?), 
"& G T 2 , has the form 

/ O kk A(#) \ 
A(0) = _ , (2.8) 

for some k X (v — k) matrix A($). Here Qjfc is the zero j x k matrix. We define the matrix 

u=[ )=u-\ 

where I k is the identity k x k matrix. By a direct calculation, one can verify that 

UA(d)U- 1 = -A(tf), 
which yields that a(A(i))) is symmetric with respect to 0. 
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Conversely, since A(0) is the Laplacian on T and its spectrum is symmetric with respect 
to 0, r is bipartite due to the main property 4) of the Laplacian. ■ 

In the following theorem we show unitary equivalence of Laplacians on graphs with multiple 
indices. 

Theorem 2.2. Let Tq = (Vq,£q) and T' = (Vq,£q) be fundamental graphs of periodic graphs 
T and T' , respectively, satisfying the following conditions: 

1) T and T' are isomorphic, i.e., there exists a bijection ip : Vq — > V ' that preserves the 
adjacency of vertices; 

2) there exists n6N such that for any e = (u, v) G Ao 

r (e)=nr(e / ), where e' = (<f(u), (p(v)) G A' . (2.9) 
Then the Laplace operators on V and V are unitarily equivalent. 

Proof. Denote the Laplace operators on T and T' by A and A', respectively, and the Floquet 
matrices for T and V by A (■#) and A' (■#) , $ G T 2 , respectively. Since r and T' are isomorphic 
and the indices of their edges satisfy the identity (12.91) . 

A(tf) = A'(n#), Vtf G T 2 . 

Then using (11.61) we obtain 

UAU- 1 = — / A(tf)dtf = — / A'H#. 
(2tt) 2 J T 2 (2tt) 2 J T 2 

If we make the change of variables = ra?, G T /2 = M 2 /(27mZ) 2 = [—7m; 7m] 2 , then we 
rewrite the last identity in the form 

UAU- 1 = — — - / A'(rv&)d-& = 7 ^ / A'ftfW = — — - / A'(tfW = U'A'U'- 1 

(27r) 2 > / T2 1 ; (27m) 2 J T 2 1 ; (27r) 2 7 T2 1 ; 

for some unitary operator U'. Thus, A and A' are unitarily equivalent. ■ 

3. Spectral analysis of Laplacians 
Below we need the following representation of the Floquet matrix A($), d G T 2 : 

A(tf) = A + A>), A = 7-^ / A(0)tft?. (3.1) 

(2tt) 2 J T 2 

From (13. ip . ( II. 7p it follows that the entries of the matrix A($) = {Ajk($)} have the form 



A ifc 



(0) = — — J2 e*< T W'*>. (3.2) 

V J e=(« jl D fc )e-4 



r(e)#0 



Proof of Theorem 11.21 i)-ii) Due to ( 13. ip we have A($) = A + A(i?), where the matrix 
A(i?), G T 2 , has the form 

A(^) = diag(A 11 ,...,A w )(^). 
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The identity ( 13. 2 j) for A 3 j has the form 

= - E cos < r ( e )' *>> 

x,- z — ' 

r(e)#0 

which yields 

_ A,,(^ ) < A^(tf) < A„(0), Vj G N„. 
Then A(tf ) < A(tf) < A(0) and we have 

A(0 O ) = A + A> ) < A(0) = A + A(0) < A + A(0) = A(0). (3.3) 
Then Proposition 16 . Xl ii gives 

An(^o) < \n(#) < A n (0), V (0, 7l) G T 2 X N„. (3.4) 
Thus, A+ = maxA n (i9) = A n (0), A~ = minA n ($) = A„($q). The last statement of the item i) 

follows from the definition of flat bands. 

iii) Since F is bipartite, the spectrum of the Laplacian on T is symmetric with respect to 
zero. From item i) it follows that Ai(0) ^ . . . ^ A^(0) are the upper endpoints of the spectral 
bands. Then — Ai(0) ^ . . . ^ — X u (0) are the lower endpoints of the spectral bands. Thus, the 
endpoints of the spectral bands A„, n G N u , are the eigenvalues of the matrix ±A(0). ■ 

Lemma 3.1. Let all edges (vj, Vj.) G Aq, 1 ^ j, k ^ v — 1, of the fundamental graph T have 
zero indices. Then 

i) The Floquet matrix A($) has the form 

A W = ( y'U ail] ) • < 3 ' 5 > 

where for each $ G T 2 the entry y{$) G C^ 1 is a vector and a(d) is a real number, A is a 
self-adjoint {y — 1) x {y — 1) matrix not depending on 

ii) If A has an eigenvalue \i with multiplicity 2, then fi is a flat band of the Laplacian A 
on the periodic graph T. 

Proof, i) This follows from Theorem ll.ll iv and the self-adjointness of A($). 

ii) Let the matrix A have a multiple eigenvalue \i. Due to Proposition 16. X I iii. there exists an 
eigenvalue A(i?) of A($) satisfying fj, ^ A(i?) ^ \i for all i? G T 2 , which yields A(-) = n — const, 
i.e., jj, is a flat band of the Laplacian A on V. m 



Proposition 3.2. Let T be obtained from the fundamental graph So of the square lattice S 
by adding u—1 (u ^ 2) vertices t>i, . . . , tv-i and u — 1 unoriented edges (vi, v u ), . . . , (i^-i, Uj,) 
with zero indices (see Figure^), v u is a single vertex of So. Then the spectrum of the Laplacian 
on r has the form 

<r(A) = a ac (A) U a fb (A), a fb (A) = {0}, (3.6) 

where the point is a flat band of multiplicity v — 2 and the absolutely continuous part cr ac (A) 
has only two spectral bands o~\ and o~i given by 

o- ac (A) = a 2 U a u *i = -<7 2 = gl,l]. (3.7) 
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"1 
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cr 2 H (7i 
I 1 — • — I 1 



Figure 5. a) The periodic graph T; b) the fundamental graph To; only two 
loops in the vertex v u are bridges; c) the spectrum of the Laplacian. 



Proof. The fundamental graph To consists of v ^ 2 vertices t>i,t>2, • • • } v u ; v — 1 unoriented 
edges (ui,i>„), . . . , (tv-i,-^) with zero indices and 2 unoriented loops in the vertex v u with 
the indices (±1,0), (0,±1). All bridges of r are loops and the graph T is bipartite. Then, 
by Theorem ll.2l iii. the spectrum of the Laplacian is completely defined by the eigenvalues of 
A(0). According to (12. 4p we have 



A(0) = 

Using the formula (16.11) we obtain 
det (A(0) - XI v ) 



( 


, 1 

\ y^+3 









1 



/u+3 
1 

i 



\ 



/u+3 



3 / 



v \v—2 



u+ 3 



■((z/ + 3)A 2 -4A + l-z/). 



Then the eigenvalues of A(0) have the form 

Ai(0) = l, A 2 (0) = . . . = A„_i(0) 



A,(0) 



z/ + 3 



Thus, due to Theorem II .2[ iii. the spectrum of the Laplacian on T has the form (13. 6p . (13. 7\i . 



V 5 V X V 2 



(b) 



a 3 fx 2 cr 2 fii o x 



-1/2 



1/2 



Figure 6. a) Graph T obtained by adding two vertices on all edges of the 
square lattice S; b) the spectrum of the Laplacian. 
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Proposition 3.3. Let T be the graph obtained from the square lattice S by adding N vertices 
on each edge ofS (for N = 2 see Figure^. Then the fundamental graph ofT has v = 2N + 1 
vertices and the spectrum of the Laplacian on T has the form 

a(A) = a ac (A)Ua fb (A), (3.8) 

where the absolutely continuous part o ac (A) = [—1, 1] consists ofN + 1 non- degenerate spectral 
bands o~i, . . . , Ojv+i, 

r 7m ~) 
( t /6 (A) = {cos ¥TT : n = l,...,iV), (3.9) 

where each flat band is simple. There are no other flat bands. 

Proof. The case N — 1 will be considered in Corollary 15.41 

Let N ^ 2. The fundamental graph To has v = 2N + 1 vertices. The matrix A(?9) is given 
by (I3.5p . where the 2iV x 2N matrix A and the vector y(i}) have the form 

/ 1 ... \ 
1 1 ... 



i J A N ®NN \ a ' 

1 :; ~ ! n> NN a n J ' AiV - 2 



1 



/ 



(3.10) 



= ( yi (A)M^)) , V.&.) = ^(l,0 l ...,e i *-) , * = 1,2, 



a(i?) = 0. According to (16. 4p all eigenvalues of A ^ have the form 

[i n = cos — - - , n = l,...,iV, (3.11) 

and they are different. Then the matrix A has N different eigenvalues fi n of the multiplicity 
2. Thus, due to Lemma I3.11 ii. the Laplacian on T has at least N = [~] flat bands fi n ,n e Nat. 
We describe cr ac (A). Identity ( 16. ip yields 

det (A(0) - A/,) = det (A - A/ 2JV ) ( — A — — XIw)-^)). (3.12) 

By a direct calculation we get 

det (A- \I 2N ) =V 2 N (X), V N (X) = det(A N -XI N ), (3.13) 



(A - A^at)- 1 = I ^ °™ ), 5 = (A* - XI N )-' , (3.14) 



y*(<&)(A - XI 2N y y($) = yl(A)B yi (A) + y* 2 (^)By 2 ^ 2 ), (3.15) 

If ( 1\ N - 1 \ 

yl(^ByM = w ^[V N _ 1 (X)+(--) cost?. J, s = l,2. (3.16) 

Substituting (13TT3]) . (I3TT5]I and (jBTTB) into ( 13TT2"]) we obtain 

det (A(#) - AJ„) =P JV (A)(-AP JV (A)--P JV _ 1 (A)- (--) A(tf)J, (3.17) 

where 

= COS^i + COS7?2- 

From the identity (16.31) it follows that 

V N (X) = -XV N ^(X) -±V N . 2 {\). 
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Using this formula we rewrite ( I3.17P in the form 

det (A(tf) - A/„) =V N (\)(V N+1 (\) - -V N ^(\) - (--j A(0)). (3.18) 
The identity ( 16. 2 p gives 

Vn ^--^ ' (3 - 19) 



-A = 2 cosv?- 



where ip is given by 

The formulas (I3.18P and ( I3.19P give 

det (A W - U,) = V N (X) ( 1) " +1 ( 3i " (jV + 2V - ^ - A W ) 
v 7 \ 2 / \ sm <^ sin if J 



= V N (X)i-j (2cos(JV + l)^-A(7?)). (3.20) 

Then the eigenvalues of the matrix A($) are defined by 

£V(A) = 0, 2cos(iV + l)ip = A(ti). 

The first identity gives all flat bands of A defined by (13. lip . From the second identity and 
the fact that the range of the function is [—2, 2] it follows that for any AG [—1,1] there 
exists "d G T 2 such that 

2cos(iV + l)cp = A($), (3.21) 
i.e., each A G [—1, 1] is an eigenvalue of A($) for some ■d G T 2 . Thus, cr ac (A) = [—1, 1]. ■ 

In order to prove Theorem ll.3[ we need the following lemma. 

Lemma 3.4. Let V = {Vjk} be a self-adjoint v x v matrix such that \Vjk\ < oo. Then 
the following estimate holds true: 



(3.22) 



- B < V < B, 

V 

where B = diag{I?i, . . . , B v }, Bj = \V jk \. 

fc=i 

Proof. Let a — (j, k), j, k — 1, . . . , v. We have the identity 

V= J2 H*)Va= E G(o), Q{a)= l -(l{a)V a + P{a)V a ) 1 (3 . 23) 

a=(j,k) a=(j,k) 

where the matrix 1(a) = {I mn (a}}, a = (j, k) is given by 



1, if (m, n) = a 
0, if (m, n) 7^ a 
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For any a = (j, k) the following estimate holds true: 

Q(a) = l(l(a)V a + P(a)V a ) ^\Q(a)\, 

W I (3 - 24) 

where \Q{a)\ = y/Q(a) Q*(a) = + I{k, k)). 

Summing (I3.24p . we obtain 



1 V V 

V= Q(<>)^ E l^ fl )l= E ^\V a \{l(j,j) + I(k,k))=J2l(j,j)J2\ V ^ = B ^ 

a=(j,k) a=(j,k) a=(j,k) j=l k=l 

which yields V ^ 5. Moreover, this yields —5 ^ V, since —V ^ B. m 



Proof of Theorem II .31 i) For each -d G T 2 the estimate (I3.22p yields 

-5(tf) ^ A(tf) <tB{ti), 5(tf) =diag{5 1 (tf),...,5 I/ (tf)}, (3.25) 

V 

B i <$) = Y^0)\- (3-26) 



k=l 



For all (j, k, #) G N 2 x T 2 it follows from (££2) that 



|A J - fc (i?)|<A ifc (0)) = — (3.27) 
The estimate (I3.27P yields that the entries Bj{d) defined by (I3.26P satisfy 

and then 

B(d) < 5(0), Vtf e T 2 . 
Using this estimate we rewrite (I3.25P in the form 

V 

- 5(0) ^ A(#) ^ 5(0), 5(0) = diag{B 1 (0), . . . , 5,(0)}, 5,(0) = ^ A jfe (0). (3.28) 

fe=i 

We use some arguments from |Kul0j . Combining (13. ip and f !3.28j> . we obtain 

Aq-B(O) ^ A{f) sC A + 5(0). (3.29) 



Thus, Proposition I6.1l ii implies 

A„(A - 5(0)) < A; ^ A n (#) < A+ s= A n (A + 5(0)), V e T 2 (3.30) 

and then 

V V 

|<r(A)| < - A;) < J2 ( A "( A ° + " A "( A o " = 2Tr5(0). (3.31) 

n=l rt=l 

Identity in (13.271) and relations f !3.28p and (13.31 j) give 

V V V j 

|a(A)| *C 2Tr5(0) = 2^5,(0) = 2 £ A, fe (0) = 2 £ , 

j=l j,k=l j,k=l v J 

which yields fll,16p . In the proof of ii) we show that the estimates (I1.16P are sharp. 
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ii) Consider the graph T shown in Figure Proposition 13.21 gives 

On the other hand, we estimate |c(A)| using (jl.lfip . For the graph T we have 

Xl = . . . = x v _i = 1, k v = v + 3. 

The fundamental graph T has only 4 oriented bridges, which are the loops in the vertex v u . 
Thus, 

£w = 4; b jk = 0, V(j,k)etft\(u,u). 
Then the estimate H 1 . 1 6 [) for the graph T has the form 

|a(A)K-^. (3.33) 

Thus, (I3.32p and (13.331) show that for the graph T the estimate (11.16P becomes an identity. 
The last statement of the item is a direct consequence of (I3.32p . 

iii) For each G T 2 the matrix A(#) has v eigenvalues, where v is odd. Theorem ll.ll vii 
gives that the spectrum er(A(i?)) is symmetric with respect to 0. Then G cr(A(0)) for any 
•d G T 2 . Therefore, \i = is a flat band of A. ■ 

Proof of Theorem 11.41 i) This follows from the facts that the point 1 is never a flat band 
(see Theorem ll.ll v) and 1 G er(A). The number of open spectral bands of the operator is 
v — r. Some of them may overlap. Then the number of gaps between them is at most u — r — 1. 
ii) This statement is a direct consequence of Propositions 13.21 13.31 15.51 ■ 

4. Perturbations of square lattice 

4.1. Fundamental graphs with one vertex. Let a fundamental graph r = (Vq, Sq) consist 
of one vertex v and any number of edges. We note that in this case the index r(e) of each 
edge e G A coincides with the vector of its coordinates in the basis a%, a 2 (the periods of T) 
and all edges of r are loops. We consider the spectrum of the Laplacian A and describe 
all isospectral graphs (with one vertex in the fundamental graph) on which the spectrum 

(7(A) =<7 0C (A) = [-1,1]. 

Since v — 1, we deduce that A($) is a scalar function given by 

A(0) = — c °sMe), G T 2 , (4.1) 



X 



v 



eeA 

where r(e) = m(e), T 2 (e)) is the vector of the coordinates of e in the basis ai,a2- The 
spectrum of the operator A on the graph T has the form 

(x(A) = a ac (A) = [A-,1], A- = minA(0) < 1. (4.2) 

Theorem 4.1. Let the fundamental graph T = (Vq, £q) consist of one vertex v of the degree 
x„ . Then 

i) er(A) = [— 1 ? 1] <^ the graph V is bipartite. 

ii) er(A) = [—1, 1] if one of the following conditions holds true: 

a) 7~i(e) is odd for all e G Aq; 

b) r 2 (e) is odd for all e G Aoj 
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c) 7~i(e) + r 2 (e) is odd for all e e Aq. 

Proof, i) The spectrum of the Laplacian consists of one spectral band. The point — 1 G c(A) 
iff the graph is bipartite (see the main property 4) of the Laplacian). It gives the required 
statement. 

Item ii) can be proved using Theorem 11.21 and the formula (11. 15)) . But we give another proof 
(by contradiction). A graph is bipartite iff there are no cycles of odd length in it (see p. 105 
in [Ur62j). Let one of the conditions a) - c) hold true and let the graph V be non-bipartite. 
Then there is a cycle with edges ei, . . . , ej G A of odd length J in it. It gives the identities 

r a (ei) + ... + r a (ej) = 0, s = l,2. 

But it contradicts all conditions a) - c), because the sum of an odd number J of odd terms is 
not 0. Thus, the graph is bipartite and item i) gives cr(A) = [—1,1]. ■ 

We consider the spectrum of the Laplacian on the square lattice S = (V,£), where the 
vertex set and the edge set are given by 

V = Z 2 , £ = {(p,p + ei ),(p,p + e 2 ) \/peZ 2 }, 

the orthonormal basis ei,e 2 coincides with the periods ai,a 2 of S, see Figure [7]a. The fun- 
damental graph So of the the square lattice S consists of one vertex v and two unoriented 
edges-loops e x = e 2 = (v, v), see Figure [7^. It is known that the spectrum of the Laplacian A 
on S has the form a (A) = cr ac (A) = [—1, 1]. 




Figure 7. a) The square lattice S; b) the fundamental graph So- 

We describe the spectrum under the perturbation of the graph S by adding some edges 
to its fundamental graph. 

Theorem 4.2. (Perturbations of the square lattice). Let S' = (!?,£') be a perturbed 
graph obtained from the square lattice S by adding N unoriented edges to its fundamental 
graph S and let A' be the set of all oriented edges of the fundamental graph ofS'. Then the 
spectrum cr(A') of the Laplacian A' on the perturbed graph S' satisfies: 

i) cr(A') = [—1,1] Ti(e) + T 2 (e) is odd for all e e A' , where (ri(e), r 2 (e)) is the 
vector of the coordinates of the edge e. 

ii) Let we add one oriented edge, i.e., two oriented edges e, e and let r = (r 1; r 2 ) be the 
vector of the coordinates of the edge e, where T\ + r 2 is even. Then cr(A') = [A _ (r), 1], where 
A~(t) satisfies 

- 1 < A~(r) < — cos < - 1 . (4.3) 

ml + \ T 2\ + 1 2 
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Moreover, the following asymptotics holds true: 

6 Irl 2 It 



71-2 Oi 1 ) , , 

\- ( T ) = -I + —— + -±^2. as |r|-»oo. (4.4) 



In particular, ifn = t%, then 

A~(r) = - min (2cos<y9 + cos2riy9). (4.5) 

3 ^e[o,7r] 

Proof, i) Let the edge set A' Q consists of edges e 1 ,e 2 ,e 1 ,e 2 of the fundamental graph of the 
square lattice and N additional unoriented edges, i.e., 2N oriented edges each of which also 
has an odd sum of the coordinates. Then by Theorem 14. H ii (the condition c), er(A') = [—1, 1]. 

In order to prove the converse we will use the proof by contradiction. Let 7~i(e) + T 2 (e) be 
even for some edge e 6 A' . Then there exists the cycle 

ei, . . . ,e ll e 2 , ■ . . , e 2 , e 

ri(e) times T2(e) times 

of odd length in the graph S'. Thus, the graph is non-bipartite and Theorem I4.1l i yields 
a(A')^[-l,l]. 

ii) Without loss of generality we may assume that ^ n ^ r 2 . From (14. ip and the fact that 
x v = 6 we deduce that 

A'(tf) = ~(cOS0i + COS ^ 2 + COS^! + T 2 $ 2 )). (4.6) 

O 

We will show (14. 3p . Using ( 14. 6 p we have 

A~(t) = min A'(i?) ^ min A'(i?) = - min (2cos<y9 + cos(t! + T 2 )ip) ^ cos^o, (4-7) 

l9eT 2 i?£T2 : 3 V56[-7T,7r] 7 

where <£> is a solution of the equation 

cos<£>o = cos(ti + T 2 )(f , (f G M. (4.8) 
The solutions of this equation have the form 

27m 

y? = ■ — , neZ. 

Ti + r 2 ± 1 

We take the solution ip given by 

7r(Tl + T 2 ) 7T 
Tl + T 2 + 1 Tl + T 2 + 1 

i.e, the nearest to n. Then this identity and (14. 7p give 

7T 1 

A (r) ^ cos^o = -cos ■ — - ^ -- . 

Ti + T 2 + I 2 

Thus, the inequality (14. 3p has been proved. 
If Ti = t 2 = 0, then ( 14. 6 p gives 

A~(0) = minA / (0) = --. (4.9) 
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Let now T\ = r 2 7^ 0. Differentiating A'($) given by ( 14. 6 p at n = r 2 we obtain the necessary 
conditions for a minimum of A': 

sin $2 

sin i?i = sin d 2 ■ 

Using the second condition in (14.101) we obtain two cases. 
Firstly, if $1 = $2, then the function A'($) has the form 

A'(tf) = ^(2cos^i + cos2ri^i). (4.11) 

Secondly, if + = 7T, then A'Od) = 1 . But inequality (14.31) gives that min A'(i?) ^ — | . 
Therefore the minimum point of the function A'(i?) is on the line $1 = $2 and thus, 

A~(t) = min A'($) = - min (2 cos-$i + cos 2ri?9i ). 

i9£T 2 , iS»i=i? 2 3 t?i6[0,7r] 

This and (jUD yield flOj) . 

We determine the asymptotics (14. 4p . We have t 2 — > 00 as |r| — > 00, since ^ tj ^ r 2 . Note 
that A'(-tf) = A'(#). Hence 

A"(r) = minA'(tf) =minA'(0), 9 = {-tt < X < vr, i?i < tf 2 ^ tt}. (4.12) 

^eT 2 tfee 

We introduce the local coordinates e (see Figure [8]) by 

7? = (tt, tt) - e, e = (ei, e 2 ) e = {0 < £l < 2tt, ^ e 2 ^ e 1 }. 
Then, using that t\ + r 2 is even, from (14. 6 j) we obtain 

A'(tf) = A'(0,r) = -( - cos ex- cos e 2 + cos /3(e)), (4.13) 

where 

/3(e) = nei + r 2 e 2 . 

Firstly, we show that the function A'($, r) achieves its global minimum on the torus T 2 in 
the domain (see Figure [S]) 

Q = {e= {s l ,£ 2 ) G [0, vr] 2 : < /3(e) ^ tt ^ tt}. (4.14) 

The function A'($, r) at the minimum point (£1,62) satisfies 

-(smei - n sin /3(e)) = 0, 



(9ei 3 
dA 1 



( sin e 2 - r 2 sin /3(e)) = 0, 



that is equivalent to 

sin Ei sin e 2 



de 2 3 

sin/3(e). (4.15) 



n t 2 

From the fact that ^ T\ ^ t 2 and the first identity in (14. 15|) it follows that sin E\ and sin e2 
at the minimum point (ei,£2) have the same sign. Thus, the minimum point (ei,ea) G [0,tt] 2 . 
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Figure 8. The local coordinates (£1,62) belong to Cl; the shaded domain fl 
contains a global minimum point of the function A'(i?,t). The point A = 

(4,4)- 

We show that in the domain Q there exists a global minimum point of the function A'(i?, r). 
Let A'(i?, r) achieve its global minimum at a point e = (§1,82) G [0,7r] 2 \ (see Figure [8]). 
Show that there exists a point e° G such that 

A'(£°,t) ^ A'(e,r). (4.16) 

We denote 

= /3(e), /3° = arccos(cos (5) G [0,vr], (4.17) 

where the branch of arccos is fixed by the condition arccos(O) = |. Note that /3 > 7T, since 
e G [0, tt] 2 \ Q. Thus, 

g = jG[0,l). (4.18) 

We define the point e° (the point A in Figure |HJ) by 

e°=(e 1 ,e 2 )=qe. (4.19) 

Using (HT7|) - (14719]) . we have 

< e° = ge s < £ s < 7r, s = l,2, 

= =/3°, (4.20) 

and hence e° G f2. Furthermore, the identity ( I4.13P gives 

A'(e°, r) = -( — cose? — cos e 2 + cos f3(e )) — -( — cos(gei) — cos(ge r 2 ) + cos/3°) 

^ - ( — cos £1 — cos e 2 + cos = - ( — cos §1 — cos £ 2 + cos ) = A' (e, r) . 



Here we have used the second identity in (14.171) . (I4.20p and the following simple inequality 

cos x ^ cos(gx), V (x, q) G [0, it] x [0, 1], 

since cos is monotonic on the segment [0, tt]. Thus, (I4.16P holds true, i.e., the function A'($, r) 
also achieves its global minimum at the point e° G Q. Then we rewrite (14.121) as 

\~(t) = minA'(#,T). (4.21) 
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Secondly, we show that at the global minimum point e = (ei, e 2 ) G of the function A'($, r) 

( £l) £ 2 ) = as \t\ -> oo. (4.22) 
|r| 

The minimum point £6(1, hence 

/3(e) ^tt, 0^e s , s = l,2. 
Since ^ n ^ r 2 , the last inequalities give 

7T 

t 2 £ 2 ^ Ti£i + r 2 e 2 = /3(e) ^ 7r, ^ e 2 < — , 

T 2 

which yields the asymptotics for the second component in (14.221) . From this asymptotics and 
the first identity in (I4.15P it follows that 

Tl . 0(1) 

smei = — sme 2 



r 2 \t\ 

and hence E\ — > or E\ — > n as |r| — > oo. But in the second case 

1 (— 1) T1 

A~(t) = — ( — cosei — cose 2 + cos(r 1 e 1 + r 2 £ 2 )) — > cosr 2 e 2 as |r| — > oo. 

3 3 

Since this contradicts the estimate (14. 3p . we conclude that e\ — > as |r| — > oo and the 
asymptotics for the first component in (14.221) holds true. Thus, (I4.22p has been proved. 

Thirdly, we obtain the asymptotics for (3(e) at the global minimum point e as \t\ — > oo. 
The second identity in (I4.15P and the asymptotics for the second component in (I4.22p give 

sin/3(e) = ^, (4.23) 
\ r \ 

which yields 

(3(e) = T\E\ + r 2 £ 2 = 7m(r) + . . , for some n(r) G Z. 

M 

Since the minimum point £6(1, n(r) = or n(r) = 1. If n(r) = 0, then 

A~(t) = — ( — cos^i — cose 2 + cos(ri£i + r 2 £ 2 )) — > — - as |r| — >■ oo. 

o o 

This again contradicts the estimate (14. 3p . Thus, the function A'(i9, r) achieves its global 
minimum on the curve 

/3(e) = nei + r 2 £ 2 = vr + ^ . (4.24) 

l T l 

This curve comes arbitrarily close to the line Ti^i + r 2 £ 2 = 7r when |r| is rather large. 

Finally, we obtain the asymptotics for the value of A'($, r) at the minimum point e = (si, e 2 ) 
as |r| — > oo. The identity (I4.24p yields 

COs(Ti£i + T 2 E 2 ) = COS ^7T + J^Pj = ~ l + ■ ( 4 - 25 ) 

Using the asymptotics (I4.22I) and (I4.25P we rewrite the identity (I4.13P as 

v ' ; 3\ 2 2 t / 6 r 4 v ; 
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This function achieves the global minimum at the point e = (ei, 62) of the curve ( 14.24)) . such 
that the square of the distance e\ + e\ from the point (0, 0) to e is minimal. Therefore, e\ + e\ 



is equal to the distance from the point (0, 0) to the curve ( I4.24p . The distance d from a point 
(x , y ) € M 2 to a line Ax + By + C = 0, A, B, C G R, is given by 



\Ax + By + C\ 

VA 2 + B 2 V ; 



Using this formula we obtain 



+ + (428) 



Substituting fl428|) into (OB|l . we have 

A- (r) = mm A'(^, r) = -1 + + ^ 
Thus, the asymptotics f)4.4p has been proved. ■ 




Figure 9. Graph S', the edges of its fundamental graph S' are marked by bold lines. 

Remark. 1) Item i) describes all iso-spectral perturbations. 

2) In ii) we estimate the end point A~(r) of the spectrum of the operator on the perturbed 
square lattice. In the case when the added edges has equal coordinates we determine A - (t). 

3) The asymptotics (14. 4p shows that if we add a "very long" edge to the fundamental graph 
of the square lattice, then the spectrum of the Laplacian is almost unchanged. 

4) For the graph S' (on Figure [9]a) cr(A') = cr ac (A') = [—1, 1], since S' is obtained from the 
square lattice S by adding the loop with the coordinates (2, 1) (and its inverse loop) to the 
fundamental graph S , 2 + 1 = 3 is odd. 

The graph S', shown on Figured, is obtained from S by adding the loop with the coordinates 
(1, 1) to So- Then Theorem I4.21 ii yields 

1 1 

A"(r) = - min (2cos<^ + cos2<^) = -- and a(A') = a ac (A') = [-1/2, 1]. 

3 ^e[o,7r] 2 
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5. Perturbations of Hexagonal lattice 

5.1. Fundamental graphs with two vertices. We consider a bipartite periodic graph T. 
By Lemma [6.4l ii. there exists a bipartite fundamental graph T = (Vq,Sq). Assume that T 
consists of two vertices V\, v 2 and any number of edges. In this case the vertices V\ and v 2 has 
the same degree x and there are no loops in T . 
For each d G T 2 the 2x2 matrix A(i?) defined by (11.71) is given by 

A (')=(a i2 I)")- = ^ E (5.1) 

e=(vi,v 2 )eAo 

The diagonal entries of A(i?) are zeroes, since there are no loops on rV For each G T 2 the 
eigenvalues Ai , A 2 (i9) of the matrix A(i9) have the form 

Ai(t?) = -A 2 (i?) = |A 12 (i?)|. 

Recall that the spectrum cr(A) of the Laplace operator on a bipartite graph with two vertices 
in the fundamental graph is symmetric with respect to and consists of two spectral bands. 
Moreover, the point 1 G cr(A). Thus we obtain 



i) a (A) = cr ac (A) = [-1, -Ao] U [A , 1], where A = min I A i2 



tfeT 2 

2 



(5.2) 



ii) cr(A) = [-1,1] A 12 (-& ) = 0, for some $ e T 

The following statement gives the method for constructing of periodic graphs with Ao > 0. 

Proposition 5.1. Let the bipartite fundamental graph Tq consist of two vertices v±,v 2 and 
A 2 multiple oriented edges (vi,v 2 ) with indices running over all values in the set 

V = {{r lj ,T 2j )eZ 2 :j = l,...,N} 

and their inverse edges. Then 

i) The function A i2 defined by ( 15. 1\) has the form 

Ai 2 W = ^Pi(e^)P 2 (e^), 

where 

N 

P s (z) = Y,z Tsj , 1*1 = 1, 8 = 1,2. 

j=l 

ii) Let P s (z) 7^ for each (z,s) G {z G C : \z\ = 1} x {1,2}. Then the spectrum cr(A) 
is symmetric with respect to and consists of exactly two spectral bands separated by a gap 
(—An, An), where A = min I A 12 (-#) I > 0. 

Proof, i) The vertices v\ and v 2 has the same degree x = A 2 . For the function Ai 2 defined 
by (15.11) we have 



A 12 (tf) = — ^2 e i{T{e) ^ = ^ e i{T1 ^ 1+T2,d2) 



N N 

— J2e iTl ^Y, eiT2k * 2 = Jf 2 Pi^ 1 )P^ 2 )- 

3=1 k=l 
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ii) Since the functions Pi,p2 have no zeroes, A 12 ($) = ^ Pi(e v ^ 1 )p2(e 1 ^ 2 ) ^ for any 
■d G T 2 . Due to (15.21) . the spectrum er(A) is symmetric with respect to and consists of 
exactly two spectral bands separated by a gap (— A , A ), where A = min Ai 2 (i?)| > 0. ■ 




FIGURE 10. A graph T; only edges of the fundamental graph Tq are shown. 

Example of a graph such that a (A) consists of two spectral bands separated by a 
gap. The functions 

P x {z) = P 2 {z) = z- 1 + l + z 2 

have no zeroes on the unit circle. Let the bipartite fundamental graph r consist of two 
vertices Vi, v 2 and 3 2 = 9 oriented edges (v±, t> 2 ) with the indices 

(0,0), (-1,0), (2,0), (0,-1), (-1,-1), (2,-1), (0,2), (-1,2), (2,2) 

and their inverse edges (see Figure [TTJ|) . Then by Proposition 15.11 the function A 12 defined by 
(15. ip has the form 

A 12 (tf) = ^ P 1 (e^)P 2 (e i ^) = (1 + e-** + e 2 ^) (l + + e 2 ^ 2 ), 

and the spectrum cr(A) on the graph T is symmetric with respect to and consists of exactly 
two spectral bands separated by a gap (— A , A ), where A > 0. Direct calculations yield 

A = min|A 12 (tf)| = | A 12 (tf*, 0*) | « 0,04; 
where #* G [0, tt] is defined by cos$* = , s = 1, 2. 

5.2. Hexagonal lattice. We discuss the Laplacian on the hexagonal lattice. 
Proposition 5.2. The Laplacian on the hexagonal lattice G satisfies 

A(0) = H D (t) + 0(\t\ 2 ) as \t\->0, (5.3) 

t = (t 1 ,t 2 )eR 2 , ^ = -1(^ + ^2), i2 = ^^_ 1?1 + 1?2 + ^ j (5.4) 
where H D {t) is the 2D Dirac operator given by 

H D (t) = oiti + <7 2 t 2 , °"i = ( ? n ) ' °" 2 = 
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Proof. The Floquet matrix A($) for the hexagonal lattice has the form 



A 12 (tf) 
It is easy to show that 



A(0)=( ), A 12 (7?) = i(l + e i * 1 +e i *»), 0=(0i,0 2 ). 

a 12 (^) o y 3 



A U (0) = O ^ tf = ±tf°, ^° = «,^)= (f ) GT 2 . 
The Taylor expansion for the entry A 12 ($) about the point $° is given by 

A 12 (0) = I(l + e **i + e **») = i (1 + e i,? S (1 + ^ - ^0) + e ^ (1 + ^ 2 - ^°)) + 0(|^ - ^°| 2 ) . (5.5) 
Using the identity e ±l ir = |(— 1 ± z\/3), we rewrite (I5.5P in the form 

tf°| 2 ) = h-it 2 +0( 

(5.6) 



2 V 

A i2 (i?) = l+l [(-l+2v / 3)(l+^i-O-(l+2V / 3)(l+^2-^ 2 )l+0(|«°| 2 ) =ti-it 2 +0(|t| 2 
3 6 



Thus, 



A C)=l (l+ ° it2 (l o it! )+o(| t | 



which yields (I5.3p . Finally, we note that the Taylor expansion for A 12 ($) about the point — 19° 
is given by the same asymptotics (15.61) . but t 2 is defined by t 2 = ~0i + 02 — ■ ■ 

The result similar to the asymptotics (15. 3p for the 2D hexagonal lattice was described 
earlier by Wallace |W47j and Slonczewski- Weiss [SW58]. Namely, they predicted that such 
a 2D monolayer material should present the branching points of the electron spectrum (Ap- 
points), where a degeneracy of the valence and conductivity band states takes place, the 
electronic state dispersion law asymptotically has a form of the double-napped cone, while 
the quasi- wave vector approaches the if-point. The equation set for the electron states in the 
vicinity of the if -point is mathematically similar to the Dirac equation for a zero-mass particle. 
These features stem from the specific symmetry of the hexagon lattice, which has two Bravais 
sublattices, and all atoms of it are situated in identical positions. These properties of this 
2D material do not depend on approximate procedures used in calculation of the electronic 
spectrum. The mentioned above similarity to the Dirac equation was used by many authors 
in the study of the graphene electronic properties: bound and resonance states, electron 
scattering, conductivity and other transport coefficients, see [FKP09] . |FK10j and references 
therein. 

5.3. Perturbed Hexagonal lattice. In order to prove Theorem ll.5l we consider the function 
T given by 

= jl + e^ 1 + e^ 2 | 2 , i?=(^i5 2 )6T 2 . (5.7) 

We have the following 

min F0&) = = 0, where = (2vr/3, -2vr/3). (5.8) 



We rewrite J 7 in the form 



J-(v) = 1 + 8 cos cos — cos — 

v ' 2 2 2 
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which yields 

maxJ 7 ^) = J^(0) = 9. (5.9) 

Proof of Theorem 11.51 Item i) is a direct consequence of ii) and iii). 

We show ii). Recall that the fundamental graph Go of the hexagonal lattice G consists of 
two vertices v%, v 2 , three multiple oriented edges ei = e 2 = e 3 = (vi,v 2 ) (Figure©) with the 
indices r(ei) = (0,0), r(e 2 ) = (1,0), r(e 3 ) = (0, 1) and their inverse edges. 

The graph G' remains bipartite iff the adding edge e connects the vertices v± and v 2 . Let 
we add the edge e = (v\, v 2 ) with an index r — (n, r 2 ) G Z 2 (and its inverse edge) to Go- For 
each d = ^) 6 T 2 the matrix A(i?) defined by (11.71) has the form 




A 12 (tf) 



1 



A(tf)=( _ I, A 12 (tf) = 7(1 + ^ +e i,?2 +e*< T »*>). 







4 



We define the point 



(0, 7r), if r 2 is odd 
$0 = { ( 7r ?0), if Tj is odd . 
(71", 7r), otherwise 

One can verify by a direct calculation that A^^o) — 0. From ii) in (15. 2p . it follows that 
cr(A') = [—1,1]. Thus, we have proved that if G' is bipartite then cr(A') = [—1,1]. The 
converse follows from item iii) of this theorem. 

iii) The graph G' is non-bipartite iff we add a loop with an index r = (r 1; r 2 ) G Z 2 (and its 
inverse loop) to G . Without loss of generality we may assume that we add a loop e = (v 2 , v 2 ). 
For the graph G' the matrix A(i?) has the form 

A 12 (tf) 



A(tf) 
where 



Ai 2 (t?) A22W 



A 12 (0) = -^=(l + e vdl +e^ 2 ), A 22 (i?) = |cos(r J i?), 1? G T 2 . 
V15 5 



The eigenvalues of the matrix A($) are given by 



\ s {d) = i cos(r, 7?) + (-l) s +y 1 cos 2 (r, i?) + 1 JF(tf) , s = 1, 2, (5.10) 

where is defined by (15.71) . Thus, the spectrum of the Laplacian on the non-bipartite 

perturbed graphene has the form 

o-(A) = [AT, At] U [Ar, At], Ar = minA s ($), At = maxA s (tf), s = l,2. 

tfeT 2 tfeT 2 

Using flBTTU]) and (J5J3]), O, we have 

At = maxAi(tf) = Ai(0) = 1, A 2 = minA 2 (^) ^ A 2 (0) = -~ . 

tfeT 2 i?eT 2 5 

The last inequality and non-bipartition of the graph G' yield (ll.20p . 
We will show ffL18|l . (05)1 . From f lBTTOj) and (RTBjl . <K^ . it follows that 

Ai(tf)>0, A 2 (tf)^0, V^GT 2 , (5.11) 
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and 



A s (tf°) = ~cos (y(n-r 2 ; 



-l) s+1 /2tt . 

cos [-j [Ti - r 2 



s = l,2, 



where = (2a , -2£). If n - t 2 G 3Z, then (CTTIjl . (15TT21 give 



A7 = minAi(i?) < A^tf ) = -, ^ A 2 

i9eT 2 5 



maxA 2 ($) > A 2 (tf°) = 0. 



< A7 = minAi(^) ^ Ai(tf°) = 0, A; 



(5.12) 



If n — r 2 ^ 3Z, then 

maxA 2 (tf) ^ A 2 (^°) = -i . 

This yields flTTgJ, (05)1 . 

The estimates (OTj) follow directly from pg| - (TOD)) . ■ 
Example of a non-bipartite perturbed graphene. Let the graph G' be obtained from 
the hexagonal lattice G by adding the loop e = (t> 2 , t> 2 ) with the index r(e) = (ti, t 2 ) = (1,0) 
(and its inverse loop) to the fundamental graph Go (Figure ITT]) . 




Figure 11. The non-bipartite perturbed graphene G'. 

Since T\ — r 2 ^ 3Z, by Theorem 11.51 we have 

a(A / )=a ac (A') = [A 2 ,A 2 + ]U[0,l]. 

We obtain numerically that 



A 2 « -0,627; 



A^ 



-0,106. 



5.4. Bipartite fundamental graphs with three vertices. Consider a bipartite graph T. 
It has a bipartite fundamental graph T = (Vo,£o) with parts V\ and V 2 . Let r consist 
of three vertices vi, t> 2 , V3 with degrees xi, x 2 , x 3 , respectively, and any number of edges. 
Without loss of generality we may assume that v\ G V\ and t> 2 ,t>3 G V 2 . 
For each 1? G T 2 the matrix A($) defined by (I1.7P is the following 3x3 matrix 

/ A 12 (0) A 13 (0) 
A(i?) = A 12 (0) I , (5.13) 

V A13W 

where 



A ls (tf) 



1 



=i(T(e),i»> 



2,3. 



(5.14) 



e=(«i, u 3 )e^o 

In the following proposition we determine the spectrum of the operator A on the graph V 
and formulate the necessary and sufficient conditions when cr(A) = [—1,1]. 
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Proposition 5.3. Let the bipartite fundamental graph T consist of three vertices and some 
number of edges. Then 

a(A) = a ac (A)Ua fb , a ac {A) = [-1, -A„] U [A , 1], a fb = {0}, 

A = min V|Ai3W| 2 + |A 12 W|2; (5.15) 

(7(A) = [-1, 1] ^ A 12 (0 O ) = A 13 (0 O ) = /or some tf e T 2 , 
where the functions Ai 2 , Ai 3 are defined by \5.1$ - 

Proof. For each ■& e T 2 the eigenvalues Ai(i?), A 2 ('t9) , A 3 ($) of the matrix A(i?) given by (I5.13P 
have the form 

Ai(i?) = -A 2 (tf) = v/|A 13 (i?)| 2 + |A 12 (i?)| 2 , A 3 (-) = /i = 0, 
which yields (15.15p . ■ 

Corollary 5.4. Lei T fre the graph obtained from the square lattice S by adding one vertex on 
each edge ofS. Then the spectrum of the Laplacian on T has the form 

o-(A) = a ac (A) U o>(A), a ac {A) = [-1, 1], tr /6 (A) = {0}. 

Proof. The graph T is a bipartite periodic graph with three vertices on the fundamental 
graph. Then by Proposition 15.31 the spectrum has the form 

°"(A) = o" ac (A) U o~fb, a ac (A) = [-1,-Ao] U [A ,l], a fb = {0}, 

where A = min a/| A 13 (i9) I 2 + lA^f^)! 2 . The functions A 12 , A 13 defined by (15.141) have the 
form 

A ls (^) = ^=(l + e -^). 

Due to identities A 12 (^ ) = A 13 (i? ) = 0, $ = (tt, 7r), Proposition 15.31 gives A = and 
cr ac (A) = [-1,1]. ■ 

5.5. The Kagome lattice. As an example of non-bipartite graph with three vertices in the 
fundamental graph consider the graph T shown in Figure [T2b. This graph is called the Kagome 
lattice. It is a lattice structure found in many natural minerals' molecular arrangements. 

The fundamental graph of the Kagome lattice consists of three vertices v%, t> 2 , f 3 each of 
which has the degree 4, six oriented edges 

ei = e 2 = (f i, f 2 ), e 3 = e 4 = (v u v 3 ), e 5 = e 6 = (v 2 , v 3 ) 

and their inverse edges. The indices of the fundamental graph edges in the coordinate system 
with the origin O (Figure H~2b) are given by 

r( ei ) = (0,0), r(e 2 ) = (-l,0), r(e 3 ) = (0,0), 

r(e 4 ) = (0, -1), r(e 5 ) = (0, 0), r(e 6 ) = (1, -1). 

Proposition 5.5. The spectrum of the Laplace operator on the Kagome lattice has the form 
a(A) = a ac (A) U a fb (A), a ac (A) = [-1/2; 1], a fb (A) = {-1/2}. (5.16) 
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Figure 12. (a) The Kagome lattice, the edges of the fundamental graph are 
marked by bold lines; (b) the support of the eigenfunction; c) the spectrum of 
the Laplacian. 



Proof. For each •& 6 T 2 the matrix A($) has the form 




The eigenvalues of this matrix are given by 

A 1 (^) = ^(l + V / ^)), = A 3 (tf) =//=--, V^GT 2 , 

where J 7 ^) is defined by (15.71) . From (15. 8p . (15.91) it follows that the ranges of the functions 
Ai($) and A 2 ('#) are [|; 1] and [— |; |], respectively. Thus, the spectrum of the operator has 
the form (15.161) and p = — | is a flat band. ■ 

Remark. 1) The eigenfunctions corresponding to the eigenvalue \i\ = —1/2 have the finite 
support shown in Figure [T2"b . The values of the eigenfunction in the vertices of the support 
are pointed out. 

2) As sets the spectra of the Laplacians on the Kagome lattice and on the triangular lattice 
(Figure [9]&) are the same. 



6. Appendix 

6.1. Properties of matrices. We recall some well-known properties of matrices (see e.g., 
[HJ85J and [HC96]), which will be used below. Let p(A) be the spectral radius of A. 

Proposition 6.1. i) Let A = {Ajk} and B = {Bjj.} be v x v matrices. If \Ajk\ ?S Bjk for all 
j, k e N u , then p(A) ^ p(B) (see Theorem 8.1.18 in [HJ85] ). 

ii) Let A,B be v x v self-adjoint matrices and let B ^ 0. Then the eigenvalues X n (A) ^ 
X n (A + B) for all n G N v (see Corollary 4.3.3 in [HJ85] ). 

Hi) Let B be the self-adjoint (v + 1) x {y + 1) matrix given by 

B={\ y 
\y* a 

for some self-adjoint v x v matrix A, some real number a and some vector y £ C u . 
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Let the eigenvalues of A and B be denoted by {fj,j} and {Xj}, respectively, and assume that 
they have been arranged in decreasing order 



fl u ^ . . . < pLi 



X U +1 ^ Xy ^ . . . ^ Al 



Then 



X u+ i ^ n v < \v ^ . . . < A 2 ^ Hi < Al, 

(see Theorem 4.3.8 in [HJ85] ). 

iv) Let M be a v x v matrix having the form 



M 



A B 
C D 



for some square matrices A, D and some matrices B, C . Then 

det M = det A ■ det (D - CA~ l B) 

(see pp. 21-22 in [HJ85] ). 
v) Let M v be a v x v finite Jacobi matrix given by 

( b 1 ... \ 

1 h 

M v = 



1 b 1 
16 



V 



/ 



where b = 2 cos if G (—2, 2) . Then 



det M u = (-1)' 



sin(z/ + l)(p 
simp 



which satisfies 

det M n+ i = b det M n - det M n _ 1; VnGN, 

with the initial conditions det Mo = 1, det Mi = b (see pp. 1511-1512 in [HC96J,). 
Moreover, the eigenvalues of the matrix M v have the form 

Tin 



A« = b + 2 cos 



u + 1 



V neff, 



(6.1) 



(6.2) 
(6.3) 

(6.4) 



6.2. The basic properties of periodic graphs and fundamental graphs. In this sub- 
section we discuss main properties of the periodic graphs. 

Proposition 6.2. The fundamental graph of a 1? -periodic graph is finite. 

Proof. Recall that we identify the vertices of the fundamental graph To with the vertices of 
the periodic graph from the bounded domain [0, l) 2 . So the number of vertices of T Q is finite 
(see item 1 of the definition of a Z 2 -periodic graph). The number of edges of T Q is also finite, 
because T contains a finite number of vertices, each of which has a finite degree (see item 2 
of the definition). ■ 

Below we need the following properties of an edge index. 
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Proposition 6.3. Let (u,v) G A and letp,q G 1? ■ Then the following identities hold true. 

i) t(u, v) = —t(v, u). 

ii) Let {v + p, v + q) G A. Then the index t{v + p,v + q) = q — p. 

Hi) The edge {u + p,v + p) G A and its index t{u + p,v + p) = r(u, v) . 
iv) Let r^(e) be the index of an edge e = (u,v) G Aq in the coordinate system with the 
origin 0\. Then 

r ( 1 )(e)= r (e) + [w-fe]-[w-6], where b = o6 x . (6.5) 
Proof. Using the definition ( II. 3p of an edge index, we have 

i) t(u,v) = [v] - [u] = -([u] - H) = -t(v,u), 

ii) t{v + p,v + q) = [v] + q — [v] — p = q — p. 

iii) The periodicity of the graph yields that the edge (u + p, v + p) G A and 

t(u + p,v + p) = [v] + p — [u] — p = [v] — [u] = t(u, v). 



v + q (v + q — b) 




Figure 13. The edge e = (u + p, v + q) in the coordinate system with the 
origin 0\ e = [u + p — b, v + q — b) in the coordinate system with the origin 0\\ 
ai, a 2 are the periods of the graph. 

iv) Let e = (u, v) be an oriented edge of the fundamental graph with an index r(e). Then 
by the definition of the fundamental graph and the formulas (II. 3p . (jl.4j) there is an edge 
e = (u + p, v + q) G A, where p, q G Z 2 are some integer vectors such that 

r(e) = q — p = r(e) (6.6) 

(Figure fT3l) . Recall that we identify the vertices of the fundamental graph with the vertices 
of the periodic graph from the set [0, l) 2 in the coordinate system with the origin O. In the 
coordinate system with the origin 0\ the edge e has the form 

e = (u + p — b, v + q — b). 

Using (Oil . (fLD and flfQ>D we obtain 

r (!)(e) = r (1) (e) = q + [v - b] -p - [u - b] = r(e) + [v - b] - [u-b]. 

Thus, the identity (16.51) has been proved. ■ 

Remark. 1) From item iv) it follows that an edge index, generally speaking, depends on 
the choice of the coordinate origin. Indeed, let the oriented edge e = (u, v) in the coordinate 
system with the origin O and the same edge e = (ui,Vi) in the system with the origin 0\ 
(Figured!]). Then the index of the edge e = (u, v) in the first system is defined by 

r(e) = [«]-[«] = (0,0)-(0,0) = (0,0) 
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Figure 14. The edge e = (u,v) in the coordinate system with the origin O; 
e = (ui,Vx) in the coordinate system with the origin 0\. 

and the index of the same edge e = (ui, V\) in the second system has the form 

r«(e) = [«i]- [«i] = (1,1) -(0,0) = (1,1). 

Item ii) shows that the index of the edge (v + p, v + q) G A does not depend on the choice of 
the coordinate origin O. It also means that the indices of all loops on the fundamental graph 
do not depend on the choice of the point O. 

2) Under the group Z 2 action the set A of oriented edges of the graph T is divided into 
equivalence classes. Each equivalence class is an oriented edge e G Ao of the fundamental 
graph IV From item iii) it follows that all edges from one equivalence class e have the same 
index that is also the index of the fundamental graph edge e. 

Now we show that a bipartite periodic graph has a bipartite fundamental graph. Recall 
that a graph T = (V, £) is bipartite if its vertex set V can be divided into two disjoint sets V\ 
and V2 (called parts of the graph) such that every edge connects vertices from different sets. 

Lemma 6.4. i) If v belongs to some part V\ of the graph, then v + 2p G V\ for each p G Z 2 . 
ii) For a bipartite periodic graph there exists a bipartite fundamental graph. 

Proof, i) The distance between two vertices in a graph is the number of edges in the shortest 
path connecting these vertices. It is known that the distance between two vertices from one 
part of a bipartite graph is even and the distance between two vertices from the different parts 
is odd (see p. 105 in |Or62j ). 
For the vertex v + p we have two cases. 

Firstly, let v + p G V\. Then the distance d(v , v + p) between v and v + p is even. Due to 
the periodicity of the graph the shortest path connecting v + p and v + 2p is obtained by the 
translation of the shortest path between v and v+p through the vector p. Thus, the distance 
d{y + p,v + 2p) is also even and hence v + 2p G V\. 

Secondly, if v +p G V2, then d{v , v + p) is odd. The distance d{v +p, v + 2p) is also odd and 
hence v + 2p G V\. Thus, in both case v + 2p G V\. 

ii) Let Tq = (Vo, £0) De a fundamental graph of a bipartite periodic graph T. If it is non- 
bipartite, then we consider another fundamental graph T' Q = (Vq, £' ) — Y/ (2Z) 2 . We will show 
that r' is bipartite. We identify the vertices of r' with the vertices of the periodic graph V 
from the set [0, 2) 2 . Let us divide Vq into two disjoint sets V{ and in the following way: 

v G V' s <^> v G V s for each s = 1, 2. (6.7) 

We show that each edge of the fundamental graph r' connects vertices from the different sets 
V{ and V2 ■ Assume the contrary, that there exists an edge e = (w, v ) G £' Q connecting the 
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vertices u,v from one set V{ or V^. Without loss of generality suppose that u,v G V[. Then 
by the definition of the fundamental graph V there is an edge (u + p, v + q) G £ for some 
p, q G 2Z. Since u, v G V 7 /, by (16. 7ft and item i) we have u + p, v + q G V\. Thus, the edge 
{u + p, v + q) on T connects the vertices from one part of T. This contradicts the bipartition 
of T. Thus, there is no edge connecting vertices from the same part of r' and hence r' is a 
bipartite fundamental graph of T. m 
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